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We describe a model of electron transfer reactions affected by local binding to the donor or acceptor 
sites of a particle in equilibrium with the solution. The statistics of fluctuations of the donor- 
acceptor energy gap caused by binding/unbinding events are non-Gaussian, and the resulting free 
energy surfaces of electron transfer are non-parabolic. The band-width of the charge-transfer optical 
transition is predicted to pass through a maximum as a function of the concentration of binding 
particles in the solution. The model is used to rationalize recent observations of pH-dependence of 
electron transfer rates involving changes in the protonation state of the donor-acceptor complex. 



I. INTRODUCTION 

Many redox reaction in chemistry and biology involve 
bound states which are weaker than common chemical 
bonds, but stronger than intermolecular interactions in 
bulk molecular liquids. A prominent example of such 
association is binding of water molecules to solutes via 
hydrogen bonds. Since the strength of such bonds typ- 
ically varies between different electronic states of the 
solute, water association can be recorded by optical 
solvatochromismiii^i^ii^ Another example is the one of 
proton equilibria involved in biological electron trans- 
fer (ET) reactions responsible for photosynthesis and 
respiratioui^ In Photosystem II, the primary donor P^go 
is oxidized by tyrosine which changes its pK^-value from 
10 to —2 upon oxidation. As a result, it cannot hold 
onto its phenolic proton in aqueous solution releasing it 
in what is argued to be a concerted electron-proton trans- 
fer mechanismi^ Another related example is photoacidity 
when photoexcited intramolecular charge transfer lowers 
pKo for the release of a proton to the solvent 4 In view 
of the wide spread of such reactions, often referred to as 
proton-coupled ET^^i^j^ in biological systems in general 
and enzymetic reactions in particular ^lii one wonders if 
their occurrence is just a coincidence caused by ubiquity 
of acid-base equilibria in proteins, or, alternatively, pro- 
tonation/dcprotonation transitions are used by nature to 
fine-tune the ET energetics. This question, also used for 
the paper title, is the subject of this study. 

It has been argued that coupling of ET to the change 
in protonation alters the reaction free energy, often turn- 
ing uphill reactions along the ET coordinate into down- 
hill reactions along a combined electron-proton reaction 
path.— In order to make this argument more qualita- 
tive, at least one more parameter, the reorganization en- 
ergy, needs to be specified J^ii^ The picture of crossing 
parabolas used to calculate barriers of charge-transfer 
reaetionsi^ involves in fact three components: the re- 
action free energy AG representing the vertical shift 
of the minima, the curvature (2A)~^ given in terms of 
the reorganization energy A, and the horizontal shift of 
the parabolas' minima responsible for the Stokes shift, 
AA: = Xoi -Xo2 (Figure [T|). 




FIG. 1: Diagram of the free energy surfaces of ET obtained 
from three parameters, the reaction free energy AG, the free 
energy curvature (2A)~^, and the Stokes shift (separation be- 
tween the minima AX). The reaction coordinate X is the 
energy separation between the donor and acceptor electronic 
states. The activated state is given by the resonance of the 
donor and acceptor levels, X = 0, while Xoi and X02 indicate 
the positions of the minima. 



The reaction coordinate X, used to monitor the 
progress of a charge-transfer reaction, is the instanta- 
neous energy gap between the donor and acceptor elec- 
tronic energy levels The vertical axis in Figure prefers 
to free energies Gi{X) as functions of energy X. They 
cross at zero energy gap X = allowing electronic tun- 
neling (i = 1,2 refer to reactants and products, respec- 
tively). In fact, the Marcus- Hush description of ET— im- 
poses one more additional constraint requiring the Stokes 
shift to be identically equal to 2A, 



aa: = 2A f3{iSX)^) 



(1) 



where P — l/{k^T). Therefore, only two parameters, 
AG and A, are required to determine the activation bar- 
rier. The strong link between the Stokes shift and the 
reorganization energy is a consequence of the Gaussian 
statistics of X. Once the statistics become non-Gaussian, 
the free energy surfaces Gi{X) are non-parabolic and eq 
[1] does not holdii^ Here, we propose a simple model for 
a charge-transfer reactions triggering the change in the 
state of binding of some particle dissolved in the solvent. 
In this development, we are less concerned with the de- 
tails of a particular binding mechanisms, but more inter- 
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FIG. 2: Diagram of particle B exchanging between the bound 
and dissolved states with the rate constant kb- The equi- 
librium population of the bound state assumes that the un- 
binding rate is much faster than the rate of ET (A;et kt); 
Ae = £2 — ei denotes the change in the binding energy caused 
by the electronic transition. 

ested in addressing a more general question of whether 
the involvement of bound states can potentially provide a 
new mechanism of tuning the energetics of ET not incor- 
porated in the picture of crossing parabolas. We find that 
the statistics of fluctuations of the donor-acceptor energy 
gap. induced by local binding/unbinding events, are non- 
Gaussian thus violating the link between the Stokes shift 
and reorganization energy given by cq[Tl The main con- 
sequence of this result is more flexibility, compared to the 
standard picture, in tuning the energetics of ET. 

II. MODEL 

We will consider a somewhat simplified situation in 
which the bound state is fully thermalized on the time- 
scale of ET i.e. on the time tet = ^et required for the 
system to climb the top of the potential barrier from the 
equilibrium bottom of the free energy surface (/cet is the 
ET rate constant). This assumption implies that popu- 
lations of bound and unbound states follow Boltzmann 
statistics along the ET reaction coordinate, a situation 
similar to the treatment of intramolecular vibrations in 
Sumi-Marcus modelJ^ This condition can be achieved 
when the rate of release of the binding particle B, fcf,, is 
much higher that the rate of ET 

kET <C h (2) 

When this requirement is satisfied, one can use statisti- 
cal mechanics to construct a one-dimensional free-energy 
surface for the ET reaction in which the exchange be- 
tween localized and delocalized states of the binding par- 
ticle produce fluctuations of the donor-acceptor energy 
gap, just as any other solvent mode interacting with the 
transferred clcctroni^ 

The release of one particle from the bound localized 
state requires the Gibbs energy gi balancing the binding 
energy with an entropy gain from moving the particle 
into the bulk. The binding energy e,; is generally different 
in the initial, D-A {i = 1), and final, D+-A^ (i = 2), 



configurations of the donor-acceptor complex (Fig. [2]). 
The difference in free energies needed to release particle 
B to the solution is equal to the difference in binding 
energies, Ag ~ Ae = £2 ^ ei- Therefore, this quantity 
can be estimated from the difference of the equilibrium 
dissociation constants pKa ■* : 

Ae = fceTlniOApK^ (3) 

The requirement of non-adiabaticity of ET, which we 
implicitly assume here, imposes an additional constraint 
on the rate of bound particle release. The change of the 
donor and acceptor energy levels caused by the motion 
of the bound particle should not break the Landau-Zencr 
non-adiabaticity condition close to donor-acceptor energy 
resonance: 

-\VnA\Ve^l (4) 

Here, Vjja is the donor-acceptor electronic coupling. As- 
suming the rate of binding energy change as e ~ k^Ae 
and VjjA of the order of 1 cm~^, eq [5] requires the rate 
of release of B to be faster than 1 ns^^. This thresh- 
old rate is comparable to the rate of proton release from 
bound protonation states of protcinsi^i while the time of 
water exchange between the protein surface and the bulk 
is faster, in the range of tens of picoseconds.^ 

Once particle B is released from its bound state, it 
becomes a part of the bulk, which is the aqueous solu- 
tion or a solvent mixture for water binding or the ionic 
atmosphere for ionic (proton) association. The interac- 
tion of the bulk with electronic states of the donor and 
acceptor is relatively well understood. For water bind- 
ing, fluctuations of the dipolar polarization produce the 
thermal noise. For electrolytes, the Debye-Hiickel elec- 
tric field both shifts the donor-acceptor energy gap and 
creates its fluctuations thus producing a corresponding 
reorganization energyfS^ This ionic reorganization energy 
~ K^e^R/eg scales quadratically with the inverse Debye- 
Hiickel length K and linearly with a solute dimension R. 
Since it is also inversely proportional to the static dielec- 
tric constant of the solvent eg , the effect of fluctuations of 
the ionic atmosphere in polar solvents is commonly small 
relative to the reorganization energy of polarization fluc- 
tuations. We will therefore assume that once particle B 
is released from its bound state, it does not interact any 
more with the transferred electron and becomes a part 
of the many-particle electrostatic potential. 

This approximation allows us to apply the tight- 
binding approximation and to use the following expres- 
sion for the Hamiltonian of the donor-acceptor system in 
a polar solvent 

H,=Eo.,-{E'^\u){u\ + E'i\b){b\)*P 

1 , (5) 

^e,\b){b\ + -P^X-'*P 

Here, Eqi incorporates the vacuum energy and the free 
energy of solvation by the electronic degrees of freedom 
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of the solvent and E,"'^ denotes the vacuum electric field 
of the donor-acceptor complex in two ET states (i = 1, 2) 
and two binding states (u,b). Also, |6)(6| and \u){u\ in cq 
Odescribe the population operators for the bound (b) and 
unbound (u) states of particle B with the binding energy 
ei. One can expect only a small change in the electric 
field of the donor-acceptor complex for weak binding of 
a neutral molecule, ~ E". However, the dependence 
of the electric field on the binding state needs to be in- 
corporated for equilibria of charged particles. 

In application to protonation/deprotonation equilib- 
ria, the Hamiltonian considered here (eq [5]) is differ- 
ent from the ones used by Cukier— and Soudackov and 
Hammes-Schifferj^ In their case, the proton was consid- 
ered to move between two localized states within the 
donor-acceptor complex, while in the problem consid- 
ered here the proton is released to the bulk and loses 
any connection to electronic transitions within the donor- 
acceptor pair except for the influence of the Debye-Hiickcl 
field it becomes a part of. The different physics of the 
problem considered here demands the different Hamilto- 
nian. 

The electric field interacts with the (nuclear) dipolar 
polarization of the solvent P, and the asterisk in eq [5] im- 
plies both the tensor contraction and the volume integra- 
tion. The statistics of P are Gaussian with the response 
function x such that the solvent reorganization energy is 



A,, 



1 



-AE„,,,*x* AE^.f,, AE 



u.b 



(6) 



Although Xu.b in eq[B]are formally the reorganization en- 
ergies of the dipolar polarization field of the solvent, we 
will attach a more general meaning to them as the reor- 
ganization energies of the nuclear degrees of freedom of 
the bulk. The corresponding formal definition is easy to 
obtain from eq[5]by replacing the Gaussian polarization 
noise with a sum of statistically independent Gaussian 
fields. 

The free energies of ET are obtained from the Hamil- 
tonian in eqOby projecting all possible nuclear motions 
on the reaction coordinate X = AH = H2 — Hi : 



-l3Gi(X} _ 



= J VPTt [6{X - AH)e''^"' ] 



(7) 



Here, VP implies functional integration over the field P 
and "Tr" refers to the sum over bound, |6), and unbound, 
states of particle B with statistical weights f^^u'- 



Tr 



[...]= J2{a\f,...\a}. 



(8) 



a — b,u 



The ratio of the statistical weights gives the entropy of 
releasing particle B to the bulk, fee ln(/„//t,). 

Taking the integral and trace in eq [7] results in the 
following equation for the free energy surfaces of ET 



{X - AE^) 



-r'ln 1 



4A„ 

^10PK«- 



pB 



MX) 



(9) 



In eq[51 G" is the free energy of the donor- acceptor com- 
plex in the unbound state. The solvation free energy 
entering G" is thus of electrostatic origin and does not 
include the free energy of binding particle B. Correspond- 
ingly, the vertical energy gap AE^ in the second sum- 
mand is given as A„ + AGu for the forward transition 
1^2 and as A^ — AGt, for the backward transition 
2 — > 1. Binding of particle B is expressed by the term 
under the logarithm where function fi{X) is the ratio 
of Boltzmann factors for activation through bound and 
unbound states: 



fi{X) = exp 



-0 



{X + Ae- AE\ 



b\2 



jX AEff 
4A„ 



(10) 

Here, AE\ are the vertical ET gaps in the bound state 
of particle B which are given as A^ -f AGb {i = 1) and 
\b~AGb (i = 2). 

In writing eqlHlwe have also represented the Boltzmann 
factor for the release of particle B in terms of pKa ^ and 
pB= — log[i?] as follows 



(11) 



In this equation gi = e,; — fcBrin(/„//b) is the free energy 
of releasing the bound particle from the molecular frag- 
ment in the donor-acceptor complex (e.g. tyrosine) and 
AGi = G" — G\ is the change in the solvation energy of 
the entire donor-acceptor complex caused by unbinding 
event. The equilibrium constants pK^*^ thus reflects the 
equilibrium of the entire complex and can be modifled 
compared to the equilibrium constants of the molecular 
fragment alone. 

The free energies Gi{X) can be used to calculate the 
first and second moments of the reaction coordinate X, 



e.g. 



where 



One gets for the average 

{X), = (1 - n,)AE: + n,{AEt - Ae) 
and for the variance 

{{SX f), = 2fcBT [XbU, + A„(l - Ui)] 

+ n,{l~ni) [AE^~AE\ + Aef 

In eqs [14] and [151 

ni=\l + lOPB-Pi^- 
is the equilibrium population of the bound site. 



(12) 



(13) 



(14) 



(15) 



(16) 




2A, -Ae 



FIG. 3: Free energy surfaces of ET (eg I17p at different val- 
ues of pB indicated in the plot assuming pKi^^ = 10 and 
pKi^^ — —2 (equilibrium constants of tyrosine). The final 
charge-transfer state (i = 2) is always deprotonated while 
the initial charge-transfer state (i — 1) is protonated at 
pB<pKi^'' and is deprotonated at pB>pKi^''. The Stokes 
shift for the protonated state 1 is 2As — Ae = 2.71 eV and 
is equal to 2 As = 2 eV for the deprotonated state 1. The 
solvent bulk component of the reaction Gibbs energy reaction 
is AG = -0.1 eV. 



Equations [THTHI provide a general description of ET 
when binding affects both the statistics of the energy 
gap fluctuations and the electronic density responsible 
for the electric field. In order to reduce the number of 
independent parameters, we first neglect the second ef- 
fect assuming that the electric field docs not change with 
binding, i.e. A„ A;, = A^, AE'f = AEf = AE,, and 
Gf = G\ = Gi. This situation most closely corresponds 
to binding of a neutral molecule (water) while the case 
of a charged particle (protonation) is postponed to the 
next section. 

The free energy surfaces of ET can then be written as 
follows: 



G,{X) = G, 



{X - AE, + Ae)2 



4A., 



(3-Hn( IQP^-'-P^-l-cxp 



(3Ae 
2Xs 



{X - AEi + Ae/2) 



(17) 



where Gi refers to the free energy of the donor acceptor 
complex with particle B released to the solution. 

Figure [3] illustrates the change of the free energy sur- 
faces with pB according to cq [17] (pKa values of tyro- 
sine). For the reaction 1^2, the system needs to climb 
the activation barrier from the bottom of the free energy 
well Gi{X) to the crossing point a.t X — 0. The bottom 
of the potential well rises with increasing pB at pB<pKi^'' 
thus resulting in a smaller activation barrier. The reac- 
tion rate depends on pB. The barrier height stops chang- 
ing once pB reaches pK^a"^ , and the reaction rate is in- 
sensitive to pB at pB>pKi^'. Notice that the rate of the 
reverse transition 2 — > 1 remains unchanged in the whole 
range of pB values. This invariance makes the effect of 
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pB 



FIG. 4: Stokes shift, AX, and the reaction coordinate varl 

\2 



ance, {{SX) )i, vs pB for a bound state with the dissociation 
constants of tyrosine (pKi^' = 10 and pKi^^ — —2) and the 
solvent reorganization energy As = 1 eV. 



pB on the ET rate distinct from the effect of the driv- 
ing force which alters the activation barriers for both the 
forward and backward reactions. 

When binding does not affect the electric field of the 
donor-acceptor complex the equations for the first and 
second cumulants (eqslTHand fTS)) simplify to 



and 



{X), = AE, - Aen, (18) 



{{5Xf), = 2k^TK + Ae'n,{l - n,) (19) 



Several important observations follow from eqs [18] and 
[T^l First, the binding/unbinding fluctuations break the 
connection (cq[T]) between the reorganization energy from 
the free energy curvature 



A, = {Pl2){{5Xf), 
and the Stokes shift 

AX = 2Xs + AeA?7., An = ^2 - ni 



(20) 



(21) 



This effect is the consequence of the local character of the 
unbinding events contrasting with the quasi- macroscopic 
nature of the bulk fluctuations resulting in eq [TJ Sec- 
ond, the reorganization energy in eqs [TO] and [20] arising 
from binding fluctuations depends on the state of the 
donor-acceptor complex, i.e. Ai ^ A2. This condition 
implies that the free energy surfaces Gi{X) must be non- 
parabolioi^ to comply with the energy conservationji^ 
G2{X) = Gi{X) + X. The effect of thermafized bound 
states on ET cannot therefore be accounted for within the 
Marcus-Hush formalism. We also note that the variance 
of the reaction coordinate in eq [TH] does not follow the 
fluctuation-dissipation theorem^! requiring {{5X)'^)i oc T 
and instead has a more complex temperature dependence 
arising from the equilibrium dissociation constant enter- 
ing the equilibrium population in eqs [15] and 1191 

The dependence of both the energy gap variance and 
the Stokes shift on pB are shown in Figure [1] The equal- 
ity of /3(((5X)^)i and AX, expected from eqs [TJ is seen 
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when both the initial and final ET states are in the same 
binding state. They are shifted by the binding energy 
Ae when the the binding states are different in the two 
ET states (left corner in Figure [4|) . The most interest- 
ing region is when pB~pKa and the term proportional to 
ni(l — rii) in cqs 1151 and 1191 can potentially dominate the 
energy gap variance. Our model thus makes a simple, ex- 
perimentally testable prediction that the band- width of a 
charge-transfer optical transition (absorption for binding 
to the donor) is expected to pass through a maximum as 
a function of pB. 



III. PROTONATION AFFECTING ELECTRON 
TRANSFER 

Measurements of oxidation rates of primary pair PjigQ 
by tyrosine in Photosystem II have produced intrigu- 
ing results 3^ The reaction rate increases with increas- 
ing pH until it levels off at pH approximately equal to 
pKa of phenolic proton. Similar results were reported 
by Hammarstrom and co-workers for intramolecular ET 
from tyrosine to Ru(III) covalently connected in a donor- 
acceptor complex i^^i^^ In order to explain the observed 
pH-dependence, they used the Marcus equation for the 
activation free energy in which a pH-dependent redox po- 
tential was substituted. The pH-dependence of the rate 
arises in their analysis from the dependence of the tyro- 
sine reduction potential on the concentration of protons 
in the solution. 

This practicC )^^'^^ which had not been anticipated 
in the original Marcus formulatiouji^ was criticized by 
Krishtalik^S and, more recently, by Saveant and co- 
workers)^ Krishtalik argued that the Gibbs free energy of 
ET reactions cannot possibly depend on the ideal mixing 
entropy of protons in the bulk, which is the origin of the 
pH-dependent term in the Nernst equation for the redox 
potentiali^ We can only add to this, absolutely correct, 
argument that using a pH-dependent AG in the Marcus 
equation clearly violates the Franck-Condon principle. 
The vertical energy gap AG 4- A does not involve entropy 
change since the nuclei do not move on the time-scale 
of electronic transitions. When the entropic pH term 
appears in AG and is not compensated by a correspond- 
ing term in A, the unphysical entropy of protons mixing 
appears in the vertical transition energy. The present 
model allows us to account for the pH-dependence of the 
ET rate without introducing unphysical approximations. 

For protonation affecting ET, B=H and pB=pH. This 
problem is, however, potentially more complex than 
binding/unbinding of neutral molecules. The process 
of deprotonation proceeds by formation of a geminate 
pair (on the time-scale of picoseconds for photocxcitcd 
states^) followed by a slower diffusion of the released pro- 
ton in the Coulomb potential of the deprotonated neg- 
ative charge4i^ This slow process may imply that the 
unbound proton will not be able to sample all possible 
states available to a particle in an ideal solution on the 
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FIG. 5: Charge shift rates from tyrosine to Ru(III) from ref 
[2^ (points) and the fit of tlie data using eq[l7] (solid lines). 
The lower curve is obtained using Au = 2.3 eV, \b = 1.8 eV, 
and the non-ergodicity multiplier a = 0.27 (AGu = —0.54 
eV and AGi, = 0.2 eV are taken from ref 1281). The upper 
curve was obtained with the same activation parameters by 
increasing the rate preexponent (see text for discussion) . The 
vertical arrow indicates the pKi^^ value for tyrosine at which 
the ET rate becomes independent of pH. 

time-scale of ET, tet- A full account of such effects^^ 
presently appears hard to achieve. We will therefore in- 
troduce an empirical non-ergodicity multiplier a to ac- 
count for the effects of insufficient sampling. The term 
pKa"*— pH in eq[9]is replaced with a(pK^*-' — pH), where 
a is a non-ergodicity coefficient between zero and unity. 
With this correction, the model qualitatively reproduces 
the dependence of the ET rate on pH observed by Sjodin 
et ali^ for the oxidation of tyrosine by Ru(III). 

Sjodin et al^ have monitored the recovery of Ru(II) 
from Ru(III) produced by flash photolysis and acceler- 
ated by electron transfer from tyrosine covalently linked 
to the ruthenium complex. Their observed rates, mon- 
itoring the arrival of electrons, mathematically corre- 
spond to projecting the complex, potentially multi- 
coordinate^ii^i^ dynamics of the system onto one sin- 
gle ET reaction coordinate, which is exactly the sce- 
nario considered here. Since proton is a charged par- 
ticle, we need the full formulation given by cq [9] to an- 
alyze the experimental rates. The reaction Gibbs ener- 
gies and reorganization energies of ET will potentially 
be different for the bound and unbound ET pathways 
and indeed the reaction free energies are AG„ = —0.54 
eV and AGf, = 0.2 eV for the unbound and bound pro- 
ton states, respectively)^ The unknown parameters are 
the rate preexponent, reorganization energies A„^b, and 
the non-ergodicity multiplier a. This latter parameter 
is expected to be small because of the slow rate of pro- 
ton dissociation from tyrosino-S and thus a lower extent 
of modulation of the donor-acceptor gap by unbinding 
events than it would be possible for full thermalization. 

In order to see if the model can qualitatively account 
for the experimental observations, we have used the rate 
preexponent, reorganization energies A^, and a as free 
parameters to fit the data from ref [2^. The result is 
shown in Figure \5\ The fitting curve (A„ = 2.3 eV, 
Xb ~ 1.8 eV, and a = 0.26) of logjo(A:ET) vs pH rises lin- 



6 




FIG. 6: Stokes shift, AX, and the variance, ({SX)^)i, calcu- 
lated from egs 1141 and 1151 using the reorganization parameters 
obtained by fitting the experimental data in Figure [S] 



early with the slope a when pH is below pK^^^ and levels 
off after reaching this value. The reorganization energies 
Xu^b here include reorganization of classical intramolecu- 
lar vibrations in addition to solvent reorganization. The 
parameters extracted from the fit of the rates result in a 
dramatic violation of eq[l]as is shown in Figure [6l 

Fitting the experimental jump in the rate at pH = 
pK^^-* requires a higher value of the rate preexponent as 
is shown by the fragment of the curve obtained by using 
the same parameters as in the low-pH fit, but varying 
the preexponent. The increase of the rate at pH = pK^^' 
has been addressed by Carra et al^ The low-pH portion 
of the curve corresponds to electronic transition accom- 
panied by simultaneous proton release (proton-coupled 
ET~). The nonadiabatic matrix element in the rate pre- 
exponent then includes small Franck-Condon overlap be- 
tween proton bound and free statesi^ This overlap disap- 
pears in the flat portion of the plot at pH > pK^^-* when 
tyrosine is deprotonated in both ET states. Only elec- 
tronic coupling enters the rate preexponent in this regime 
resulting in a higher rate. A more quantitative analysis 
would require extensive calculations. Since the model 
presented here does not aim at a detailed quantitative 



analysis of proton-coupled ET, we limit our discussion to 
qualitative observations only. 



IV. CONCLUSIONS 

Traditional theories of ET in polar liquids emphasize 
activation of electronic transitions by long-range, quasi- 
macroscopic solvent fluctuations.— Local binding to the 
donor and acceptor sites provides an additional modula- 
tion of the donor-acceptor energy gap. This study ad- 
dressed the question of whether this additional thermal 
noise can be accounted for within the standard frame- 
work of Gaussian fluctuations of the energy gap, that 
is by adjusting the magnitudes of the driving force and 
reorganization energy. 

We have found that the statistics of binding events are 
non-Gaussian, and the resultant free energy surfaces can- 
not be reduced to crossing parabolas. The model predicts 
a regime of a significant dependence of the activation bar- 
rier on the concentration of the binding particles in the 
solution. When the concentration pB is close to the bind- 
ing equilibrium constant pK^, the variance of the energy 
gap passes through a maximum which can be observed by 
spectroscopy of charge-transfer bands. Despite these new 
features, the main effect of binding is in shifting the free 
energy gap of ET, as was suggested previously;^ leaving 
the reorganization energy mostly unaffected and within 
the realm of standard models. Finally, the model helps 
to rationalize some recent observations of the dependence 
of rates of intramolecular ET, involving deprotonation of 
reactants, on the pH of the solution. 
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